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ABSTRACT 

Turbulent viscosity is believed to circularize and synchronize binary orbits and to 
damp stellar oscillations. It is also believed that when the tidal period is shorter than 
the turnover time of the largest eddies, turbulent viscosity is partially suppressed. 
The degree of suppression, however, is disputed. We re-examine both of these beliefs 
via (i) direct perturbative calculations, linearizing the fluid equations on a turbulent 
background; and (ii) numerical integration of a chaotic dynamical system subject 
to periodic forcing. We find that dissipation of rapid tides is severely suppressed. 
Furthermore, circularization of late-type binaries does not occur by turbulent 
convection if it occurs on the main sequence. 
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1. Introduction 


Eddy viscosity in convection zones is presumed to dissipate the shear associated with variable 
tides, and hence gradually to circularize and synchronize close binaries [ Zahn 1966 , Zahn 1977 , 
Zahn 1992 and references therein]. This theory of tidal circularization has been tested against 
binaries containing giant stars, and the agreement is satisfactory (Verbunt fe Phinney 199G|). 


This paper will focus on a theoretical controversy concerning eddy viscosity whose outcome 
may have consequences not only for tidal circularization, but also for stellar pulsations and 
oscillations. Assuming a Kolmogorov cascade of eddies, we have on scales A smaller than the 
mixing length (l that 


v\ 


T\ 


v\ 



1/3 

vi, 

- 2/3 

n, 

2/3 


(i) 









- 2 - 


where vi, i\, and vi are the convective velocity, turnover time, and effective viscosity on the scale 
of the mixing length, which is comparabale to the pressure scale height, H; while v\, t\, and v\ 
are the corresponding quantities on smaller scales A <1. If the eddies on scales are space-filling, 
then the total turbulent viscosity vt ~ zy • But when the period (P) of the tide or oscillation 
is less than 77 , momentum transport by the large eddies is inhibited. During half a period (or 
perhaps one should consider P/2n), fluid in such an eddy travels no farther than (P/2r/)Z. Zahn 


(1966, 1989) and |Zahn Sz Bouchet (1989)| therefore suppose that the eddy viscosity is reduced by 
the same factor; that is, they take 


v\ = min (P/2t a , 1). 

Under this hypothesis, the largest eddies continue to dominate the total viscosity, 
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Goldreich &; Nicholson (1977, henceforth GN) and Goldreich fc Keely (1977) have taken a 
different view. These authors argue that the viscosity v\ on any given scale should be severely 
suppressed—perhaps by an exponential factor—when P -C 2ttt\, so that u t is dominated by the 
largest eddies whose turnover time is less than P/2ir. Since v\ oc r\ [cf. eq. (;lj, it follows that 


ZAp = — Vil min 


m 2 ,i 

\ 2 t m) 


( 4 ) 


GN’s brief argument can be quoted in its entirety: 


To appreciate this, it is necessary to recognize that turbulent eddies have lifetimes 
which are comparable to their turnover times. Thus, even though the large convective 
eddies move across distances of order /tt/t c in a tidal period, they do not exchange 
momentum with the mean flow on this timescale. The contribution to zAp made by the 
largest eddies is likely to be smaller than zap max by at least one additional factor of 


Underlying this intuition, perhaps, is an analogy with integrable hamiltonian systems, in which 
secular absorption occurs only through resonances; if such a system is perturbed rapidly compared 
to any of its autonomous frequencies, no resonance with the perturbing forces is possible. 
Convective turbulence is neither integrable nor hamiltonian. |Goldman &; Mazeh (1991) have 
motivated GN’s reduction factor by means of an analogy with the kinetic theory of gases, in which 
they equate 77 to an intermolecular collision time. Unfortunately, as shown in §4, Verbunt fc 


Phinney 1991 ’s calibration against giants cannot be used to decide between Eqs. @ and 


Another argument in favor of GN’s prescription can be made by considering the time inverse 
of absorption. Consider a convecting star in an initially exactly circular, short-period orbit. 
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Convection causes small fluctuations in the density distribution within the star and corresponding 
departures of the exterior potential from symmetry; the latter modulate the binary orbit, albeit 


at a very low level (Phinney 1992). If r; <C P, however, the energy of the orbit is adiabatically 
invariant, so that no secular transfer of energy from the eddies to the orbit can occur. It would 
therefore seem to follow that the large eddies should not absorb energy from a slightly eccentric 
short-period orbit. Arguments from time-reversibility are very powerful for thermal systems, but 
turbulent convection is not reversible because of the flow of energy from large to small scales. So 
this argument is not conclusive. 

What is at stake? 


Tidal circularization. In solar-type stars, 17 ~ 20 d in the middle of the convection zone. Thus 
2ttti is substantially longer than the orbital period of circularized main-sequence binaries. 
Furthermore, the turnover time increases with the depth of the convection zone both on the 
main sequence towards lower masses, and in the pre-main-sequence phase towards larger 


radii. Zahn &; Bouchet (1989) therefore conclude that pre-main-sequence circularization 
would be ineffective if the more severe reduction factor (^) were correct. 

Stellar pulsation The fundamental pulsation period in variable stars is (tf/GM) 1 / 2 , which 
is usually shorter than 77 unless the convection zone is thin. Gonczi (1982)| has argued that 
the red edge of the instability strip is more consistent with Zahn’s prescription (g) than with 
GN’s (g). 

Solar oscillations Goldreich and collaborators have constructed a theory for the excitation 
and damping of the solar p-modes by convection | Goldreich fc Keely (1977)1; |Goldreich fc| 


Kumar (1988) ; Goldreich, Kumar, fe Murray (1994) and references therein]. This theory 
is generally in good agreement with the observed energies of the modes, certainly more so 
than any alternative, but it would have to be severely modified if eq. (]|) were applicable. 
Acoustic emission by turbulence can be estimated, at least roughly, with some confidence; 
if the damping processes were enhanced then the energy in the convection would be lower 
than observed. 


The outline of our paper is as follows. In the next section we attack the fluid equations 
directly, deriving linearized equations for the perturbation in the turbulent velocities due to the 
tide. The in-phase correlation of the perturbed velocity with the tidal force determines the rate of 
tidal work done on the turbulence. In §3 we study a toy model for tidally perturbed convection: a 
system of strongly coupled nonlinear oscillators driven by a periodic external force. By numerical 
integrations, we find the rate of absorption of energy of this system as a function of the frequency 
of the external forcing. In §4 we discuss the implications of our results for circularization of 
late-type spectroscopic binaries. 
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2. Perturbative Methods for Eddy Viscosity 

We divide the total velocity field into a large-scale flow Vt associated with the equilibrium 
tide (or stellar pulsation) and a small-scale flow v c associated with the convection: v = Vt + v c - 
The small-scale motions are approximately incompressible and isentropic, so we describe them by 
the eulerian equations 

d t v c + V t ■ W c + v c ■ VV t + v c • Vu c + Vw c = 0, (5) 

V • v c = 0. (6) 

The quantity w c incorporates pressure and gravitational accelerations, which are gradients of 
scalars in the present approximation. Buoyancy and thermal diffusion are not represented here, 
although they are necessary to sustain convection; nor do these equations incorporate coriolis 
forces or molecular viscosity. We have considered these complications at length. The conservative 
effects (buoyancy and coriolis forces) do not seem to affect our conclusions at all, although they 
complicate the algebra, and the true dissipative terms are important only if they are large enough 
to damp the tide directly. 

We adopt a local approximation in which the lengthscales of the convection are infinitesimal 

compared to those of the tide. Locality is inherent in the notion of an effective viscosity and is 

responsible for much of the usefulness of the concept. In this approximation, the tide couples to 
the turbulence through the first-order spatial derivatives of Vt. Equation (]6|) contains the term 
V t ■ Vu c , however, in which Vt is not differentiated. This term includes bulk transport of the local 
eddies, which does not directly contribute to dissipation. We eliminate this term by adopting 
a semi-lagrangian coordinate mesh that deforms to follow the large-scale displacement. Thus if 
{x z \i G 1,2,3} are fixed coordinates, then the semilagrangian system {x /J } is defined implicitly by 

(£),-*<*•*>■ <7) 

Since we expect that only first spatial derivatives of Vt are important for dissipation, we represent 
Vt as a linear function of the coordinates, 

V = A (t) ■ x or Vt = A}(t)x J , (8) 

in which the matrix A depends upon time but not space. The symmetric and antisymmetric parts 
of A describe the local instantaneous shear and vorticity of the tidal flow. 

Equations (@)-(|§|) can be formally integrated to obtain a time-dependent linear transformation 
between the two coordinate systems: 

x = B(f) ■ x' or x l = Bjx'i. (9) 

The matrix B(t) is related to A(f) by the differential equation 


B (t) = A(t)-B(t) 


(10) 
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which has the power-series solution 

t 


B(t) = 1 + / dtiA(ti) + j dt i j dt 2 A(ti) ■ A(i 2 ) + • • • 


( 11 ) 


if x = x' in the far past. The contravariant components of the small-scale velocity field in the 
semilagrangian system are related to the original components by 


fi a —*/ t-x — X 

v = — — -v J or v = B 


dxi 

and the equations of motion become 

d t v + 2B _1 ■ A • B ■ v + v ■ VV + V'w = 0, 

V'-if = 0. 


( 12 ) 


(13) 

(14) 


(The subscript “c” will be omitted from the small-scale quantities henceforth.) As a check, we 
note that if A were purely antisymmetric, then A • u = 12 x u for an any vector u, with Q(t) the 
instantaneous angular velocity; B (t) = (B T ) -1 would be a rotation matrix; and eq. H) would 
reduce to the familiar form 

<V + 20 ! x if + if ■ VV + V'w = 0, 

with Q' = B _1 • O. 

The average rate of work done (per unit mass) on the small-scale velocities by the tide is 


1 / d- 


1 / d- 


£ = -(—v-v) = -(—v'-B T -B-v l 

2 \dt / 2 \dt. 


= (v 1 ■ B T • B • d t v') + -(if ■ B T ■ (A T + A)B • v') 


(15) 


Actually we intend a triple average: the angle brackets denote an average over space and over 
realizations of the turbulence, while the overline denotes an average over time. Using equation 
( 0 ) and spatial integration by parts, one has 


£ = —-(u' • B T ■ (A T + A)B • v') 


= ~ 2 ^ ’ ( AT + A )u). 


(16) 


This could have been derived directly from equation (|j) without use of the primed coordinates, 
but the semilagrangian system proves convenient in the next step. Henceforth let A be symmetric, 
A = A t , since there should not be any dissipation associated with the antisymmetric (rotational) 
part of the the tidal velocity gradient. 


The velocity v or v* appearing in the dissipation formula (16) is the actual velocity field in the 
presence of the tide, not the the unperturbed turbulence uq that would exist without it. When the 










- 6 - 


tidal shear is small, it is appropriate to evaluate 5v' = v' — vo by perturbation theory in ascending 
powers of A. We can then express the dissipation rate ([LG|) in terms of the statistical properties of 
UO) even though it is beyond our power to evaluate vq directly. We assume that the unperturbed 
turbulent velocities can be described as a stationary random process with correlation function 


C™(Ax, At) = (y™(x + Ax, t + At)v 7 o (x, t )). 
We do not assume that vq (x, 2) has gaussian statistics. 


(17) 


The first-order contribution to £ is obtained by replacing v 1 with vq in equation (1C). This 
contribution vanishes, as can be demonstrated with a temporal integration by parts. At second 
order, 

£2 = ~2(v 0 ■ A ■ 5v') (18) 

To evaluate this, one must calculate dv 1 to first order in the tide. Linearizing equation (13), 


dtSif^x', t ) = —2A(t) • F 0 (x', t) — vq ■ V'di? — 5i? ■ V'vo — (pressure term). 


(19) 


The pressure term has not been written out because its only role is to enforce V' • 5v* = 0, which 
can be imposed directly in Fourier space. We will express Sv 1 , vq, and A in terms of their Fourier 
transforms: 


iut 5i7{k,LU 

iuJ % (k,u>) 

1 r 


Sif (x',t) = 

v 0 (x',t) = 

A(f) = Real A(h2)e 


dui 

f d3 k Stf 

2t r J 

(27r) 3 

dio 

f d3 k jk-g> 

2t r J 

(2t r) 3 


—iflt 


2 L 


A(0)e" int + A(-ft)e' 


iQt 


( 20 ) 


Since the quantities on the left are real, their transforms satisfy 5if(—k, — u) = 5i?(k,Lu)*, 

Vq(— k, —co) = vo(k,u)*, and A(— £1) = A(f2)*. When necessary, u> and hi are considered to have an 
infinitesimal positive imaginary part so that 5v'(x',t) vanishes as t —> — 00 . Equation (HD becomes 

— iiodv (. k , u) = —Pj: ■ | A(fl) • Fo(cu — 12, k) + A(—12) • vo(iv + 12, k) + 

/' dio’ f dk! 


27r J (27r) 3 - 


vo(k 1 , u>')dif(k — — lo') + Sv 1 (F, u>')vo(k — k',u> — oj') ■ ( ik ) > , 


where 


P f = I - k~ 2 kk 


( 21 ) 

( 22 ) 


is a projection operator that maintains incompressibility, 

k ■ 5v (k, uj) = 0. 


(23) 
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The terms on the second line in equation (^l|) are small compared to those on the first line. 
Neglecting the second line, one has 5v' ~ 1 , where V c is typical velocity for the convection. 

It follows that the second line is smaller than the first by a factor rvw' V C /(QL C ) = (flr c ) 1 , where L c 
is a typical convective lengthscale, and r c = L c /V c is the corresponding eddy turnover time. Note 
that we now have two characteristic dimensionless parameters: the first is the tidal strain, fl _ 1 |A|, 
which is always small in detached binaries; the second, (f 2 r c ) -1 , is completely independent of the 
first and may or may not be small. When the latter is small, we may formally solve equation (^l|) 
as a power series in (f7r c ) _1 . Let Si :n i? denote the contribution to Si? that is first order in the 
tidal strain and n th order in (flr c ) _1 . Successive terms are generated recursively: substituting 
J2k= 1 <5i ,kv' for Si? in the second line of (21), one solves for Y?k=\ ,k'?- 

To lowest order in (17 t c ) -1 , we have 


(k , uj) = - Pr ■ A(fi) ■ vq{lo — fl,k) + A(—fl) ■ vq{uj + fl, k) 

UJ L 

and upon substituting this into equation (|l^), we have 

duj f d 3 k (i 


£2,2 = 


We have made use of 


J 2tt J (27r) 3 \uj 
C(jfe,w-n) ■ A*(fi) -Pj-A(fi)]. 


TVace[C (k,u + n) ■ A(fi) ■ P r ■ A*(Q) + 


(vo(k', uj')vq (k, uj)} = (27r) 4 h 3 (F + k)6(u> + u/)C(/c, a;), 


(24) 


(25) 


(26) 


where C(k,ui) is the Fourier transform of the correlation function (|l7|), and the transpose symbol 
is used for clarity: vv T is a matrix with components v l v? The delta functions in equation ( f2G| ) 
result from stationary statistics: that is, from the fact that the correlation function ( |l7| ) depends 
only on the differences between the two points at which the turbulent velocity is correlated but 
not on the absolute positions and times. 


The subscripts follow the scheme described above: that is, corresponding to the term hgn-i 
in the expansion of the perturbed velocity, there is through equation ( |l 8 |) a contribution £2 ,n to 
the dissipation rate that is second order in the tide and n th order in the convection. In addition 
to two powers of A and n powers of vq. £ 2,7 1 also contains (n — 2) spatial derivatives applied to vq 
and (n — 1) time integrations. Assuming that each differentiated factor of vq is of order t ~ 1 , and 
that each time integration donates a factor fl -1 , we may expect that 


4n ~ |A| 2 F c 2 r c (fir c ) 1 - n . (27) 

Since the hrst nonvanishing term is £ 2 , 2 ; the above scaling agrees with the prescription (^) of 
Zahn. It will now be shown, however, that £ 2,2 is actually much smaller than predicted by the 
scaling (|27|). 
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The trace in equation ( |25|) is real because 

Trace(C ■ A ■ Pj ■ A*) = TVace(C A A*f = Trace(A* • Pi ■ A • C*) = Trace(C* • A* • Pi • A), 


in which the arguments of the various matrices have been temporarily suppressed, and we have 
made use of A = A r , P? = Pt and C{k,uj) T = C (k,uj)*. On the other hand, the term iu 
in eq. (25) appears to be purely imaginary. Since £2 is real by definition [eqn. (IS)], it would 
seem that £ 2,2 must vanish. As noted above, however, all frequencies should be endowed with 
an infinitesimal positive imaginary part, because one wants the particular solution for which Sv 1 
vanishes as t —> — 00 . Hence one must interpret i/ui as 


lim- 

e—>0 00 + ie 


lim 


lUJ 


■*0 ^ 


+ 


lim — 7 . 

e—>0 co 2 


+ € 


2 ' 


(28) 


The first term on the right is indeed imaginary. The second term is real and represents it5(uj). 
This term gives a nonzero result for £ 2 , 2 , which can be written using A(— £1) = A*(H) as 


£2,2 - 2 


d 3 k 

(2i r) 3 


Trace[C(fe, +H) • A(fi) • P^ • A*(H) + C(k, —£l) ■ A(-H) • P : • A*(-fi)]. (29) 


We can compare the last result with the dissipation rate provided by a true molecular viscosity 
12 : 

£ visc = -v' Trace[A(H) • A*(fi)]. (30) 

In general, the effective viscosity provided by equation (|29[) is anisotropic. If the turbulence is 
isotropic—though there is no good reason that convective turbulence should be—then 


C(k,u) = C(k,Sl)Pft (31) 

with C(k, — H) = C(k, fl). In this case, equation ( |29| ) can be simplihed to 

7 r d 3 k 

£ 2,2 - —Trace[A(H) • A*(ft)] / -j^C{k, £1). (32) 

15 J (27rj J 

We have taken Trace(A) = 0, since all of our derivations assume incompressible motion; our 
arguments would need to be generalized to treat damping of compressive motions, such as solar 
p modes. The integral in equation (f32|) is the kinetic energy of the turbulence per unit frequency 
per unit mass at the frequency of the tide. Denoting this quantity by £(£i), we can express the 
turbulent effective viscosity in the isotropic case as 

u T ~ ^S(Q), (33) 

15 

where £(uj) is normalized so that the average kinetic energy per unit mass is 


\\v(x,t )\ 2 



— OO 


(34) 
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The 


symbol acknowledges that this is only the first term of an expansion in powers of (f It c 




Whether one considers the isotropic form (32) or the more general form (|29]), the leading-order 
contribution to the turbulent dissipation rate depends entirely upon the power spectral density of 
the turbulent velocities at the frequency of the tide. Imagine an artificial turbulence in which the 
eddies on all scales had correlation times r c » 12 . Then C (k, ±12) in equation (|2^) and £(±I2) 

in equation (33) would be exponentially small, and £ 2.2 oc exp[—(I2r c )]. 


In fact, even if the largest eddies are much slower than the tide, there will be smaller and less 
energetic eddies with decorrelation times t\ < IW 1 . Our perturbative expansion is invalid for these 
small eddies, which correspond to large values of the wavenumber k in eq. (^). It seems likely, 
therefore, that the expansion in powers of (fire)” 1 is at best asymptotic rather than convergent. 


There is a question about how to translate the heuristic Kolmogorov model into a prediction 
for £{u). In the usual interpretation, the turbulent power per unit frequency at frequency to is 
due to eddies of size A such that their turnover time t\ ~ co~ 1 . The standard Kolmogorov scalings 
(p]) applied to eq. (29) then lead to GN’s prescription ([|) up to a dimensionless constant of 
order unity. Another interpretation is possible, however, and leads to a prescription intermediate 
between equations (£|) and (||). The correlation function whose Fourier transform appears in the 
dissipation rate ( |29| ) is measured at fixed spatial positions. A small eddy of size A ~ is carried 
past a fixed point by the resultant velocity of all larger eddies in which it is embedded; this velocity 
is ~ vi, the velocity of the largest eddies. The apparent frequency of the small eddy measured by 
a fixed observer is kvi, which is larger than the turnover rate 7 T 1 by a factor ~ (fcZ) 1 / 3 [cf. eq. 
©]■ According to this logic and our result (|33f), the effective turbulent velocity scales with the 
tidal period as 

Y p \ 5 / 3 

{2^1) ’ 


z^T = vd min 


(35) 


Physically, however, we do not see how the mere fact that a small eddy is borne along by a larger 
one should affect its rate of absorption of energy from a spatially uniform tidal shear. Perhaps 
higher-order terms replace the eulerian power spectrum C(fc,I2) by a lagrangian equivalent that 
measures the decorrelation rate of the turbulence following the fluid, in which case one reverts to 
the earlier interpretation, and the exponent 5/3 is replaced by 2, in agreement with eq. (!)• Even 
if eq. ( |35| ) is correct, it is closer to eq. (|]) than to Zahn’s prescription (;|). In the rest of our paper, 
therefore, we adopt GN’s prescription. 


3. Chaotic Dynamical Models 

As a toy model for convection, we have considered a chain of coupled nonlinear oscillators, 
which may loosely be considered to be interacting vortices. The lagrangian of the unforced system 
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is 


N—l 


L = X \y 2 i ~ \(yi+ 1 - y^ 3 ~ 7(^+1- y*Y 


(36) 


i =0 


with fixed end-point boundary conditions yo = yjv = 0. The linearized normal modes of this 
system have zero frequency, so their interactions at finite amplitude should be analogous to strong 
turbulence. 


We have simulated this system (36) by direct numerical integration of the equations of motion 
for N = 63. We set the initial displacement of all oscillators to zero and assign velocities drawn 
independently from a gaussian distribution of unit variance, after which the velocities are rescaled 
so that the total energy 

N -1 

I „ I „ I 

(37) 


E = X \y 2 i + livw -1 h ) 3 + \iyi +1 - Vit 


i =0 

is iV — 1 at the start of all runs. 

To verify that the system is chaotic, we estimate the maximum Lyapunov exponent (cf. 
Argyris, Faust &; Haase 1 994] ) 

a v i \d(t) 

A = Inn iog 2 


t— KX) 


d{ 0 ) 


(38) 


where d is the infinitesimal distance in phase space between neighboring trajectories, and 
the logarithms to base 2 reflect the convention that information regarding the initial state 
is measured in bits. We solve the non-linear equations of motion for a reference trajectory 
y(t) = (yi(t), • • .,y N -i(t)}, 

Y(t) = F[Y(t)], (39) 

and the linearized variational equations describing the displacement of a neighboring trajectory, 


dF 

SY W “ aF 


Y(t) 


SY(t). 


The distance between neighboring trajectories is measured by 


d(t) = ||<5T(t)|| = 


N -1 1 / 2 

i— 1 


(40) 


(41) 


Although the limit t —+ oo is inaccessible, we find the increase in d(t ) to be accurately exponential 
and estimate A ~ 2.7. 

A quantity analogous to the energy per unit mass per unit frequency of convection [eq. 
is the temporal power spectrum of the oscillator velocities, 

2 . 


£{u) = 


1 


N - 1 


N—l 




i=i 


Vi(u) = J e lujt yi{t), 
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which we compute by discrete Fourier transforms and average over many runs. The spectrum 
(Fig- 0 ) has broad peaks at low frequencies and a smoothly declining tail above a characteristic 
frequency defined by 


w dyn = 


(rn 2 ) 

(ii>'ii 2 ) 


1/2 


0.465. 


(42) 


Equipartition of energy appears to hold among the spatial Fourier modes defined by 


Qk{t) = 



(43) 


Although the eigenfrequencies of these modes vanish, their characteristic frequencies at finite 
amplitude, a; fci dyn = (<?fc> 1/2 /(<?i) 1/2 , are found to be approximately proportional to wavenumber k, 
and and the largest of these frequencies are comparable with Wdyn- 

To simulate the tide, we add the force term 


Fj = eyj cos (fit) 


(44) 


to the equation of motion of a single oscillator (j = 22). The combination ecos(fit) is analogous 
to the tidal shear A (t) of §2, and we usually take e ~ 0.05tdd yn . 


To determine the secular energy absorption rate (E), the total energy E{t) defined by eq. (| 
is recorded as a time series, smoothed with a low-pass filter, and fit to a straight line. The results 
for (E) as a function of If are compared in Figure 0 with the power spectrum (42). The energy 
absorbed from the “tide” drops off sharply when 0 > cj<j yn , as expected. 


Since our oscillator system is hamiltonian, the argument from time-reversability in §1 is 
applicable to it. That the system can absorb very little energy from a high-frequency tide therefore 
should not come as a surprise. In fact, the fluctuation-dissipation theorem is applicable here 
and predicts that energy is absorbed at a rate propotional to £[oj) (cf. |Landau fe Lifshitz 1980 , 

. It would be interesting to allow for dissipation at large k and stochastic 
forcing at small k [cf. eq. (|4^) above] in the unperturbed oscillator system; this would break 
time-reversal symmetry, establish a cascade from large scales to small, and bring the model closer 
to real convective turbulence. But to achieve a good separation of timescales would require many 
more than 64 oscillators and might be a challenging computational exercise. 


Nelson &; Tremaine 1997 


4. Application to Late-Type Binary Stars 

In this paper, we have revisited the question of how efficiently an oscillatory large-scale shear 
is dissipated by turbulence when the reduced period of the shear ( P/2tt ) is shorter than the 
correlation time of the energy-bearing eddies (r c ). By a direct perturbative approach proceeding 
from the dynamical equations of incompressible turbulence (§2), we find that the dissipation rate 
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Fig. 1.— Energy absorption rate vs. 
power spectrum eq. (34) (line). 


forcing frequency (points), as compared with the velocity 


8 (u) 





- 13 - 


is reduced by a factor scaling as ( P/2-kt c ) 2 [or perhaps (P/2ttt c ) 5 R , depending on what one takes 
for the temporal correlation of the turbulence on timescales r < t c \. The quadratic scaling is also 
found in numerical simulations of a forced nonlinear oscillator chain (§3). Our results therefore 
agree substantially with the quadratic suppression ([|) rather than the linear one ([][). 

The distinction between these prescriptions is important astrophysically because it is often 
the case that P <C 2i rr c . In the remainder of this paper, we apply our results to the circularization 
of late-type binaries. Our treatment follows that of Zahn (1989, henceforth Z89), except that we 
adopt the quadratic inefficiency factor rather than the linear one. 


4.1. The long-period limit 


To establish a baseline for later discussion, we first recalculate the circularization times of 
late-type main-sequence binaries as though their periods were long compared to all eddy turnover 
times. Even at this naive level, there appears to be a discrepancy between the tidal theory and the 
observed maximum period at which main-sequence binaries are circularized. When we take into 
account the inefficiency of slow convection, in the next subsection, the discrepancy will become 
severe. 


Consider a binary consisting of two stars identical to the present sun. The circularization rate 
follows from equation (21) of Z89: 


where 


is Zahn’s “friction time,” 


-i d , A (R 

& ^-s lne = 84 *7 (« 

2 \ !/3 


tf = 


MR 


and the dimensionless parameter A is defined by 

R 


A = 205.6 


t 


f 


MR S 


jpur r s dr , 


(45) 


(46) 


(47) 


where r c is the inner radius of the convection zone. The friction time is only very weakly 
dependent on stellar mass and spectral type; for the sun, tf = 0.432 yr. The parameter A is a 
volume average of the dynamic turbulent viscosity, p^T, weighted by the square of the tidal shear; 
the radial displacement of the equilibrium tide varies approximately as r 4 . We have multiplied 
Zahn’s expression for the circularization rate by a factor of two to account for the dissipation in 
both stellar envelopes. 


Using mixing-length theory, Z89 derives 
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P'R 
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75 AirRt 
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6 c y / 3 


a 4 /3£2/3 x -2/3( l _ x f. 


(48) 





- 14 - 


Here a is the usual mixing-length parameter (the ratio l/H of mixing length to pressure scale 
height), and x = r/R is the fractional radius. The superscript on Vj!' indicates that this estimate 
of the turbulent viscosity assumes Qr c <C 1. The dimensionless parameter c < 1 specifies how 
much of the potential energy released by convection is converted to kinetic energy on the scale of 
the mixing length; after discussing the choices made by other authors, Zahn takes c = 9/32, but 
the dissipation rate is not very sensitive to c because of the cube root. The constant E describes 
the run of density with radius in a polytropic approximation to the convection zone: 

M ( 1 /? l n r2 /I r\l 3 / 2 

" e [z {—)] ’ < 49 > 

where n ~ 3/2 is the polytropic index. Comparing the final expression above with the tabulated 
model of [Bahcall, Pinsonneault, &; Wasserburg (1995) , one finds that E ~ 6.1 at the base of the 
convection zone, which occurs at 0.712Rq, and that E decreases slowly to ~ 5.4 at the outermost 
tabulated point (0.9464Rq). 

We therefore find for the circularization time 



t c - nc « 9.0 x 10 10 


where R/a has been eliminated in favor of the orbital period. As this formula shows, t c 
sensitive to the stars’ mean density p = ‘SM/AnR^. 


(50) 


is 


For a stellar population of age t, the transition between circular and eccentric orbits should 
occur at a period P c i rc such that 


^•circ (Pci 


1 




(51) 


so that mean initial eccentricities of order 0.5 will have been reduced to a few percent or less. By 
this criterion, 

-dm 


Pci 


d 


(52) 


Recall that eq. (52) is only an upper bound because it neglects the inefficiency of slow 
convection; yet observations suggest that P c i rc exceeds this bound. The halo-binary sample 
of |Latham et al. (1992a)| indicates P c irc i> 12 d for metal-poor, high-proper-motion dwarfs 
(p 0."27 yr” 1 , [Fe/H] < —1.6). [Latham et al. (1992a) conclude P c i rc « 19 d based principally 
on a binary with period 18.74 d and eccentricity 0.043 ± 0.021. This system is a single-lined 
spectroscopic binary, so that its companion might be evolved. But the sample also contains two 
cicular double-lined systems with P = 10.74 and 11.73 d. Based on colors, metallicities, and 
proper motions, all four stars are clearly on the main sequence and well below the turnoff (Laird, 
Carney, fe Latham 1988| ). [] Also, -P c i rc ~ 10 d has been estimated for main-sequence binaries in 
the old open cluster M67 (|Mathieu Sz Mazeh 1988 , Mathieu, Latham, &; Griffin 1990 ). 


Rrom the most metal-poor 16-Gyr theoretical isochrone of Bergbusch & VandenBerg (1992), which fits M92, the 
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4.2. Inefficient turbulence 


Revising the theoretical circularization time to allow for inefficient dissipation of rapidly- 
varying tides makes the conflict between theory and observation still more acute. 


The rapidity of the tide relative to the convection is characterized globally by the dimensionless 
quantity 

2 irt f ( t f \ f 10 d\ 

~99[ - j . 


i] = 


-ftide 

which we have defined by analogy with Z89’s parameter rj = 2fy/P t ide = t// 7r - Locally at each 
radius within the convection zone, however, the convective timescale is 


(53) 


Ti(r) = — = 


where l = aH is the mixing length and is the turbulent viscosity in the long-period limit 
[eq. (|48|)1 . Where r^(r) > P ti d e /(27r), we apply GN’s prescription ([|), so that the dynamical 
viscosity becomes 

, ,, T\J 

-cx~ 3 (l - xy 1 . (55) 




4 °)’ 


(54) 


M) 1—2 
pv>p = r] 
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At: Rtf 

Taking c = 9/32 as before, we use the smaller of the expressions 


and 


the integral 


for A. Notice that the reduced viscosity (55) is independent of the mixing length and of the 
density parameter E. Where eq. (|5^) applies, in the lower parts of the convection zone, the total 
dissipation rate receives equal contributions from equal logarithmic intervals in exterior mass, 

M - M r (cf. GN). 

Figure |L2| shows the tidal parameter A [eq. (|47|) 1 as a function of rf. At r/ ^ 10, A is scarcely 


reduced below its long-period limit because the r 8 factor in eq. (47) emphasizes the outer parts 
of the convection zone where the local timescale r<tj. The asymptotic scaling A oc r / -2 is 
approached only gradually; near r / = 10 2 , A oc r/ -1 ' 5 . Another salient feature of this plot is 
the relative size of A in fully convective versus solar-type stars. Fully convective stars enjoy an 
advantage of almost an order of magnitude in the long-period limit because they are less centrally 
condensed; however, this advantage almost disappears for r] 3> 10 because the fully convective star 
has a longer local timescale 17 (r) than the solar-type star at the same fractional radius. To identify 
the point corresponding to a period of 10 days, we have assumed a fully convective protostellar 
sun with T e ff 


160. 


4070 K, following Zahn fe Bouchet (1989) , and find 7 / 

Fig. shows that for solar-type stars, t circ is increased by ~ 13(10 d/P ) 1,5 with respect to 
the estimate ([5(]). The corresponding P c i rc of old solar-type stars should be « 1.3 d instead of 


mean densities of the stars in these two binaries are approximately 3.6 and 4.8 times larger than that of the sun if 
they lie on the main sequence. Eq. ( |50| ) then predicts tcirc ~ 2.7 x 10 12 and 6 x 10 12 yr, respectively. If the stars were 
giants, then they would be 5.4 and 8.2 magnitudes more luminous, would be much more distant, and would have 
absurdly large transverse velocities. 
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Fig. 2.— Tidal dissipation parameter A [eq. (47)] versus dimensionless tidal frequency rf [eq. (|53|)j. 
Solid line: solar model. Dashed line: fully convective star. Circles: two suns (solid) and two 
protosuns (open) in a 10 d orbit (see text). 
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the value given in eq. (|52|), and still less for stars whose main-sequence lifetimes are longer than 
16 Gyr. 

This figure also implies that A(?/) ~ A(0) for the giant-star binaries studied by Verbunt & 
Phinney (1996, henceforth VP). The friction time tf oc T eff 4 ^ 3 at hxed mass [eq. ([h]), but T c g- 
varies only very slowly along the giant branch: VP quote T c fj/T eff Q ~ (R/Rq)~ 0A1 . Even at the 
tip of the giant branch (R ~ 160-Rq), 2irtf ~ 5.7 yr, whereas the binary periods of interest are 
of order 10 2 ' 5 — 10 3 ' 3 d, so that r{ 6. Thus turbulent viscosity is expected to be fully efficient 
in these binaries, and the results of our paper are consistent with the good agreement VP found 
between their giant-star binaries and the long-period limit of Zahn’s theory. 


4.3. Pre-main-sequence circularization 


In summary, one cannot explain the circularization of old main-sequence binaries by turbulent 
convective viscosity. This is true even if one ignores the reduced efficiency of turbulent dissipation 
at short tidal periods. If taken into account, as our results would indicate that it must be, the 
inefficiency widens the discrepancy between the observed and computed circularization period to 
a full order of magnitude. 




Zahn V Bouchet (1989) (henceforth ZB) have suggested that circularization occurs primarily 
before the main sequence. Although the T Tauri phase lasts ^ 10 7 yr, it may dominate the 
circularization process because the stellar radius is larger during this phase than on the main 
sequence (cf Stahler & Walter 1993). Other things being equal, the circularization rate varies as 
R 8 [eq. fl45|)]. Hence a factor ~ 3 in radius may compensate for a factor ~ 10 -4 in lifetime. 


ZB have integrated the circularization rate against time for low-mass T Tauri stars following 
the evolutionary tracks of Stahler and collaborators ( Stahler 1983 , Stahler 1988 ). According to 
ZB’s results, circularization during the pre-main-sequence (PMS) phase is effective out to periods 
of 7.2 — 8.5 d for stellar masses 0.5 — 1.25 Mq. This is based on the linear prescription (||) for the 
inefficiency factor. ZB remark that the time-averaged circularization rate is 2 — 5 times smaller 
than it would have been if they had evaluated the turbulent viscosity in the long-period limit. 
We believe, however, that the quadratic prescription ([]) should be used instead of the linear one. 
For ZB’s initial protosun (M = Mq, R = 4.79 Rq, T e fj = 4070 K), the dashed curve in Fig. |4.2| 
indicates a reduction factor of 130. Although we have not repeated ZB’s evolutionary calculations, 
we can rescale their results to obtain a new prediction for P c irc by assuming that our initial 
reduction factor (130) should replace their maximum reduction factor (5). Since A oc t/~ 1-5 
the range of present interest, eq. @) implies t c i rc oc .P 3 ' 83 , and so 

/130\ _1/3 - 83 


over 


If 


circ,PMS 


V 5 J 


8 d « 3.4 d 


(56) 


Given ZB’s assumptions about the protosun, therefore, we agree with their conclusion that 
circularization by turbulent convection is much more effective in the T Tauri phase than on the 
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main sequence, although we have halved their prediction for P c i rc 


One can imagine circumstances in which pre-main-sequence circularization would be much 
more effective. The evolution of protostars is still uncertain; perhaps protostellar radii are briefly 
larger than Stahler and ZB suppose. Tidal interaction with an intrabinary or circumbinary disk 
may be important, though in some cases the interaction may increase eccentricity (Artymowicz et 


al. 1991, Artymowicz fe Lubow 1994| ). Observations of PMS spectroscopic binaries are difficult to 
interpret. Difficulties arise from the paucity of known short-period orbits, but perhaps also from 
genuine variety in the physical circumstances of PMS systems. |Mathieu (1994) concludes that the 
evidence to date indicates P c irc ~ 4 d among PMS binaries with solar-mass components. This 
agrees, perhaps fortuitously, with our revision of ZB’s theory. At any rate, this is much shorter 
than the conservative inference -P c i rc ;> 12 d from observations of old halo binaries. 


Perhaps the PMS history of the low-metallicity halo binaries was systematically different 
from that of present-day T-Tauri binaries in such a way as to enhance PMS circularization of the 
former—e.g. longer-lived or more massive disks. This could give rise to a spurious correlation 
between P c i rc and the age (or metallicity) of the binary population. Evidence against this idea is 
that solar-mass binaries in the relatively young Pleiades and Hyades clusters show intermediate 
transition periods (7.05 and 8.5 d, |Mathieu et al. 1992 and references therein). 


Pending more data for short-period binaries of all ages, it seems that gradual circularization 
does proceed on the main sequence, but that some mechanism other than turbulent convection 
must be responsible. 


We would like to acknowledge much helpful and patient advice from Bohdan Paczynski 
concerning the theoretical evolution and observed properties of stars in general and of binaries in 
particular. This work was supported by NASA under grant NAG5-2796. 
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